Time evolution of continuous-time quantum walks on dynamical percolation graphs 
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We study the time evolution of continuous-time quantum walks on randomly changing graphs. 
At certain moments edges of the graph appear or disappear with a given probability. We focus on 
the case when the time interval between subsequent changes of the graph tends to zero. We derive 
explicit formulae for the general evolution in this limit. We find that the percolation in this limit 
causes an effective time rescaling. Independently of the graph and the initial state of the walk, the 
time is rescaled by the probability of keeping and edge. Both the individual trajectories for a single 
system and average properties with a superoperator formalism are discussed. We give an analytical 
proof for our theorem and we also present results from numerical simulations of the phenomena for 
different graphs. 
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I. INTRODUCTION 

Randomness in graph theory has been an ever develop- 
ing field since the pioneering works by Erdos and Renyi 
[I]. Today, there is a refined mathematical theory for 
treating random graphs [2 and, at the same time, sta- 
tistical methods borrowed from physics continue to be a 
motor for research in random network theory [3\ Ap- 
plications of classical random graph theory often involve 
the notion of dynamics on the random graph, describing 
phenomena ranging from disease spreading [1] to perco- 
lation [3]. Random walks provide a simple, fundamental 
model for dynamics on graphs and random graphs, being 
one of the first problems in random graph theory and 
later used throughout its applications [B]. 

Quantum dynamics, and especially quantum walks on 
graphs have been a topic of much interest in the past few 
years [7 , 8 . The concept of quantum walk was first intro- 
duced by Aharonov et al. [9] as a discrete-time dynamical 
system. The continuous-time quantum walk was shown 
to be an effective tool for quantum computation in the 
pioneering paper by Farhi and Gutmann [10] . Possible 
applications in quantum information processing triggered 
thorough research on all aspects of quantum walks. 

Apart from their role as a universal quantum informa- 
tion processing primitive [TT], continuous-time quantum 
walks are a conceptually simple, but effective model for 
transport [12] . Randomization of the underlying graph 
structure is motivated by uncontrolled processes in the 
physical model, i.e. transport in disordered media [13] . 
similar to percolation in a classical system. The percola- 
tion quantum walk problem has also been considered for 
discrete time quantum walks [14] . 

In the theory of static percolation, we have a given 



graph, in which we can keep an edge with probability 
A, or remove it with probability 1 — A. During the time 
evolution, we keep the graph constant. The ensemble av- 
erage of a quantity of the system can be calculated as an 
average over the possible realizations of the percolation 
graph. In case of continuous-time quantum walks, the 
static percolation problem has already been considered 
under the name: statistical networks [Tfjl ITH] . 

For quantum walks, changing the underlying graph 
during the time evolution randomly at some rate leads to 
a different problem. In the context of percolation, this is 
called the dynamical percolation problem 17 . Dynami- 
cal percolation for discrete time quantum walks has been 
studied numerically 14 and its asymptotic behavior for 
one dimensional systems also analytically [T5] , 

In this paper we consider the problem of continuous- 
time quantum walk on a dynamical percolation graph. 
First, we define the problem of dynamical percolation for 
continuous-time quantum walks. Unlike in the discrete 
time case, there is no natural timing for changing the 
graph structure. One has to introduce a time step for 
changes. Then we concentrate on the continuous per- 
colation limit, when changes in the graph may happen 
infinitely frequently. 

The structure of the paper is the following. In Sec- 
tion [TTJ we briefly review the basic properties of the 
continuous-time quantum walk. Then in Section III we 



consider the time evolution of the walk on a percolation 
graph. We introduce a simple formula for the time evo- 
lution and we analitycally prove it. We extend our result 
for ensemble average of quantum walks on percolation 
graphs with a superoperator formalism in Section |IV| Fi- 
nally, in Section [Vj we show some numerical examples to 
demonstrate behavior for finite parameters. 
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II. CONTINUOUS-TIME QUANTUM WALK 

Quantum walks have two main types: the discrete- and 
the continuous-time case. In the discrete-time quantum 
walk the walker moves on a discrete graph, and changes 
its position in discrete time intervals [T§1 H2U] . 

The concept of continuous-time quantum walk was de- 
fined by Farhi and Gutmann [TO]. We have a G(V,E) 
undirected graph, and let us denote the vertices with 
a = 1 . . . N. Then we construct an iV-dimensional 
Hilbert space with an orthogonal {\a}} basis correspond- 
ing to the graph where a — 1 ... N and (a\b) = 5 a .b- The 
unitary time evolution reads 

U(t) = e' i6t , (1) 
where the matrix elements of the Hamiltonian are 

{—7, if I a) and |6) are neighbouring states, 
0, if j a) and |6) not neighbouring, 
fry, if|a) = |6). 

(2) 

In the definition of H the k parameter is the degree of 
the vertex, and 7 is a time- independent constant. For 
simplicity, we can choose 7 = 1, it only rescales time. 
The wave function of the system after time t reads 



|*(t)) = e -* i "|#(o)) 



(3) 



The probability, that we measure the walker at the 
state 1 6) after time t if we started the walk from the state 
I a) can be calculated as 



*Ut)= (b\e-' lHt \a) 



(4) 



III. 



QUANTUM WALK ON A PERCOLATION 
GRAPH 



In the previous sections we briefly introduced the basic 
model of percolation, and the continuous-time quantum 
walk. In this section we consider how the dynamics of the 
continuous-time quantum walk changes if we modify the 
underlying graph randomly during the walk. We focus 
on the case when the time interval between the changes 
tends to zero. 

Let us suppose that we have a graph with N edges, and 
the nodes are connected by a given starting geometry 
(e.g. ring). Then in given time steps we change the 
underlying graph randomly. In each step we can keep an 
edge of the original graph with probability A, or discard 
it with 1 — A. We evolve the system for time T with step 
size t. In this case the number of the steps is S = T/t. 

For every possible realization of the graph, we can in- 
troduce the 



unitary operator according to ([I]), where the number of 
the r realization indexes is R = 2 . 

Let us suppose that the initial state of the system is 
\tp(0)). In each step, we act on the actual state of the 
system with an U r (r) unitary operator. Introduce the 
following notation: let U Ts (r) be the time evolution op- 
erator according to the randomly generated graph in the 
sth step, in accordance with |2]). With this notation, the 
state at time T, after S steps reads 

|V(T)) = U rs (r)...U r2 (r)U ri (T)\m) 

= e~ lil "s T . . . e -iHr 2 T e -iH ri r |^ ^ _ (g) 

Since the H r operators do not commute, we can use 
the Zassenhaus product formula: 

e t(A+B) = e tA e tB TJ e t*cf '*> + 0(^+1) 5 

lim e tA *e tA * . . . e tA ^^ . . . e'"^ = e^U A > , 

71— >00 

(7) 

where the Cj A ' B \ ci Aj ^ parameters are the Zassenhaus 
exponents [Si] [55]. We are interested in the r — >• limes, 
therefore we will neglect the higher order terms in r: 



iHr 



e -^Hr p +H rt ) + (t*). (8) 



With this approximation, if t is small enough, the final 
state reads 

\m) = (< ' r%: " +0(r 2 )) |^(0)) . (9) 

Let us consider the probabilities of the realizations. 
The probability of the realization that contains N edges 
is p( N > = X N . The probability of a realization, in which 
one edge is missing is p^ -1 ) = A w_1 (l — A), and this 
case is degenerate, because the number of the graphs 
that contain N — 1 edges is TV. It can be seen that the 
probability of a configuration that contains k edges is 
p(k) _ _ \^N-k^ an( j degeneracy is (^). 

We can decompose every H r matrix in the following 
way 



Hr 



E 



Ei- 



(10) 



where £ r is the set of the k values according to the edges 
present in the rth realization, and E% is a matrix belong- 
ing to one edge in the graph. For example, if the H r 
realization contains two edges, the first and the second 
one, then H r = E\ + E 2 , and the Hamiltonian of the 
graph without percolation reads 



N 



UAt) 



-iH r r 



(5) 



(11) 



fc=i 
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Our task is now to evaluate the sum 

s 

s = E Hr > ■ 



(12) 



If S > R then some H r matrix occurs more than once 
in §. The probability that we get the H r matrix x times 
from the S steps follows a binomial distribution. If S 3> 
R then we can use the law of large numbers. We can 
estimate the number of occurrences of one H r matrix in 
the S sum as p r S + 0(VS), where p r is the probability 
of realizing the graph that belongs to H r . The accuracy 
of our estimation can be calculated from the standard 
deviation y/ Sp r (l — p r ), it leads to an 0(\/~S) error. 

If the realization contains k edges, thenp r = p^ k \ Now 
let us calculate, how many times a chosen E% matrix 
occurs in the § sum. This matrix belongs to the k edge 
in the graph. Choose j edges, where j = 1 ... N — 1, 
and put them to the graph next to the k edge. We can 
do it in ^J 1 ) ways. Therefore the number of the H r 
realizations, that contains E^ reads 



N-l 

E 

3=0 



N - 1 
3 



(13) 



A given H r matrix occurs p r S + 0(y/S) times in the § 
sum, it means, that the sum contains the matrix 



N-l 

E 

3=0 



N - 1 

3 



= SX + 0{VS) 

times. It is true for every fc, therefore 

S 



H . 



(14) 



(15) 



§ = E Hr ° = sx + 

s=l 

With this result the exponential in ^ can be written as 

s 

-ir^ H rs = —iTH (A + 0(y/r)) . (16) 



8=1 



Taking the r — > limes, we arrive at our main result 



|^(T)) = lim JJ U rk (r) |^(0)) - U{T\) |^(0)) . (17) 

r— y 

fe=l 

The effect of a finite time step r will be examined nu- 
merically in a later section. 

Finally one can ask what happens if we take the ex- 
tremal values for the A probability. Both cases are trivial. 
If A = then we have a null graph, there is no time evo- 
lution, the system stays in the initial state. If A = 1 
then we get the percolation free graph in each step, we 
get back the percolation free evolution. We can simply 
arrive at these results if we write A = or 1 into the ( 17) 
equation. 



IV. SUPEROPERATOR FORMALISM 

In this section we would like to consider the average 
properties of the randomly evolved systems that we in- 
troduced in the previous section. The standard way to 
calculate the ensemble average of a quantity of a system 
reads 



(...) = 



1 



Ei 



t=i 



(18) 



where T is the number of the possible trajectories, and t 
denotes one trajectory |T2"] . 

In each step, we can keep an edge with probability A, 
or discard it. Let us denote the length of a time inter- 
val with r. Let us consider a periodical chain with N 
nodes, and evolve the system until time T. In this case, 
the graph changes T/t times. If N > 2, then the graph 
has N edges, therefore the number of possible graph re- 
alizations is R — 2 N . It means, that for time T the 
number of the possible trajectories is T = [2 N )^ T / T \ be- 
cause we can choose from R graphs in each step. For 
a big graph with good time resolution this explicit eval- 
uation becomes very difficult. For example the number 
of the trajectories for a periodical chain with N = 15 
nodes with S = 5000 steps is over 10 22577 (we use these 
parameters later for a numerical simulation with the su- 
peroperator formalism) . 

Instead of evaluating all possible trajectories one by 
one, we can use a superoperator formalism |23H25j . in 
which we use the formerly introduced U r (r) unitary op- 
erators. In one step, we act on the density operator with 
the 



l T (p) = Y j p r U r (T)pUl{T) 



(19) 



r=l 



superoperator, where p r is the probability of the graph 
realization denoted by r as before, and we sum over all 
possible realizations. If we evaluate the system for time 
T with t steps, then we should act S = T/t times with 
the 4>t superoperator, and the final state reads 



p(t) = MM^(- . . Mp))) = 4 s Hm) 



(20) 



If we use the |5| form of U r (r), the final state can be 
written as 



rs = l \rs-i = l 



n=i 

R 



= P ri -Pr s e^ H ^ T ...e- lH '^p(0)e lH ^ T ...e lHr s T . 

(21) 



n...rs=i 
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Here we can use the same treatment in the r — > limes 
that we use in the previous section. In that section we 
demonstrated, that in the r — > limes the typical be- 
haviour of the possible trajectories is only a time rescal- 
ing, therefore if we take the average of the trajectories 
we again get the time rescaled dynamics. Therefore the 
act of the superoperator can be written in a very simple 
form: 



p(T) = lim 



4>i : \p(0)) = U(T\)mu\T\). (22) 



V. NUMERICAL SIMULATIONS 

In the previous sections we demonstrated that, if we 
have a percolation graph in which we keep an edge with 
probability A in every r time step, then in the dynamics 
of a continuous-time quantum walk on this graph, the 
percolation only causes a time rescaling if the r time 
interval tends to zero. In our calculation we use the r — > 
limes. Of course, it is only a mathematical limes, in a 
physical system we cannot get his exact limes. In this 
section we show some concrete example to show, that for 
every time T there exists a small enough r parameter, for 
which we get a very good estimation for the final state. 
This can be relevant for experimental realization of the 
phenomena. 

In our simulations we calculate the probability that 
the walker returns to the origin. It is a perfect param- 
eter to illustrate our analytical result, it characterizes 
the transport properties of the walk. Many articles use 
the return probability to investigate the properties of the 
continuous-time quantum walk [26 -29] . 

In our first example we consider the continuous-time 
quantum walk on a percolated 2D integer lattice that 
contains N = 100 nodes. We address the nodes by i,j, 
it means the ith node in the jth column. We evolve the 
system with r = 10~ 4 time steps in S — 10 5 steps. We 
start the system from the |\&(0)) = |45, 45) state, and we 
calculate the probability of being at this node. We use 
different A values, and we plot both the numerical and 
the theoretical values according to (17). The results are 
in Fig. [I] Incidentally, we note that the shape of the 
curves are very close to the appropriate Bessel function 
[3"0] if we scale them by the T — » AT transformation in the 
percolation case. These results suggest that our results 
can be extended for infinite 2D integer lattice graphs. 

In the next example we demonstrate that the equation 



( 22 1 gives a very good approximation for the time evo- 



lution with finite r values if we take an average over the 
trajectories, i.e. we use the superoperator formalism. We 
take a periodical chain with N = 15 nodes, and we evolve 
the system with t = 0.004 steps from the p(0) = |1)(1| 
initial state in S = 5000 steps. We use different values 
for the probability, and we also plot the theoretical line 
from (22). Our result are in Fig. [5] Again, we can be 



see that the starting shape of the curves are very close 
to the Jo(2t) Bessel function [31 . This fact again allows 
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FIG. 1: The numerical simulation of one trajectories in a 
10 x 10 2D lattice in 10 5 steps with different probabilities and 
the theoretical line. 
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FIG. 2: The numerical simulation of the superoperator for- 
malism in a periodical chain with TV = 15 nodes. The param- 
eters are S = 5000 and r = 0.004. 



us to suggest that our results are valid for the infinite ID 
integer lattice graph. 

In the third figure we illustrate that the precision of 
the ( 22 ) formula can be arbitrary fine if we choose small 
enough r. We consider the system for time T = 10, and 
we divide the [0, T] interval into S part, i.e. we evolve 



the system in S step until T. Therefore if S is big, then 
the r = T/S parameter is small. In our simulation we 
use a periodical chain with N = 10 nodes with different 
initial state and A probabilities. For every S values, we 
calculate the probability of being at the |1) state both 
with the numerical simulation and the theory, and we 
take the maximum difference between them. Our results 
can be seen in Fig. [3] It is nicely seen that the difference 
between the simulation and the theory is a decreasing 
function of the number of the steps. 
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FIG. 3: The maximum difference between the numerical sim- 



ulation and the ( 22 1 theoretical values of the probability of 
being at the |1) state in an N — 10 periodical chain. The 
time interval that was divided into S parts was [0, 10]. 
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The initial state of the system was /5(0) = |1)(1|. With 
these settings we can see that after a long time (after lot 
of steps), the probability that we measure the walker at 
the initial state can be described with 



Pi(t) 



1 

N 



(23) 



This asymptotic dynamics has already been considered 
for the discrete-time quantum walk [18 , and it seems 
that similar methods can be applied to the continuous- 
time case [32) . In our simulation we calculate the prob- 
ability of being at the |1) state, it is the p\^\ element of 
the density matrix. In Fig. [4] we show the numerical 
results, clearly indicating that the probability tends to 
the theoretical value of 1/5 according to (23 1. 



VI. CONCLUSION 

In this article we studied the time evolution of the 
continuous-time quantum walk on a percolation graph. 
Both unique trajectories and their average properties by 
using a superoperator formalism have been examined. 
We found that if the t time scale of changing the graph 
tends to zero, then the time evolution can be described 
with a simple analytically formula. In fact, the percola- 
tion causes only a time rescaling compared to the perco- 
lation free case. This result is independent of from the 
underlying graph and we could prove this theorem ana- 
lytically in the r — > limes. Naturally, in a real physical 
system we can expect that r is a small but finite number. 
We have presented numerical examples to demonstrate 
how the system follows our prediction in the long time 
limit if we choose the r time steps small enough. The dis- 
cussed system may serve as a simple model for transport 
in disordered media. 



FIG. 4: The probability of measure the walker at the first 
node in an N = 5 periodical chain. We evolve the system for 
S = 10 5 steps with r = 0.01 time intervals. We choose the 
probability of keeping an edge for A = 0.5. 

In our last simulation we look at what happens, if we 
evolve the system for a long time with not too small r 
time steps. We consider a periodical chain with N = 5 
nodes for S = 10 5 steps with r = 0.01 time interval, A = 
0.5 probability, and we use the superoperator formalism. 
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